Rectification of spin currents in spin chains 
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We study spin transport in nonitinerant one-dimensional quantum spin ciiains. Motivated by 
possible applications in spintronics, we consider rectification effects in both ferromagnetic and an- 
tiferromagnetic systems. We find tiiat the crucial ingredients in designing a system that displays a 
nonzero rectification current are an anisotropy in the exchange interaction of the spin chain com- 
bined with an offset magnetic field. For both ferromagnetic and antiferromagnetic systems we can 
exploit the gap in the excitation spectrum that is created by a bulk anisotropy to obtain a mea- 
surable rectification effect at realistic magnetic fields. For antiferromagnetic systems we also find 
that we can achieve a similar effect by introducing a magnetic impurity, obtained by altering two 
neighboring bonds in the spin Hamiltonian. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Ds, 75.76.+j 



I. INTRODUCTION 

Spintronics is the field of spin-based electronics, in 
which the spin degree of freedom is taken as a carrier 
of information instead of the charge, as is the case in 
traditional electronics. The field has seen a rapid de- 
velopment over the last decades.^ Among the potential 
advantages of spin-based devices (which could be op- 
erated at high switching speed) in comparison to their 
charge-based counterparts are longer coherence times of 
the spins and lower energy dissipation associated with 
performing logical operations and transporting the infor- 
mation carriers.— The latter advantage is often strongly 
reduced by the fact that most schemes for spin-based 
devices rely on spin currents that are generated as a by- 
product of associated charge currents.— 

The focus of the present work lies on transport of spin 
excitations in nonitinerant one-dimensional quantum sys- 
tems (spin chains). In such systems magnetization can be 
transported by magnons (spin waves) or spinons (domain 
walls), without any transport of charge. Since transport 
only occurs in the spin sector the energy dissipation as- 
sociated with the transport of a single spin in nonitiner- 
ant systems'* is much lower than in itinerant systems,-^ 
in which spin transport is a by-product of the trans- 
port of charge. Furthermore, nonitinerant systems have 
also been proposed for performing low-power logic op- 
erations.^ Since energy dissipation is such an important 
limitation in modern spin and electronic devices, these 
facts make nonitinerant systems a promising candidate 
for applications in spintronics 4 

Possible candidates for spin chains are, for instance, 
bulk structures of KCuFa, SR2CUO3, or SrCu02, in 
which the exchange interaction between different chains 
in the crystal is much weaker than the intrachain in- 
teractionr^^*^ There exist spin chains with rather larger 
anisotropy in the exchange interaction, for instance, 
CS2C0CI4 (see Ref. fisl ). In a recent experiment^ - it has 
been shown that the cuprate SrCu02 has a mean-free 
path for spinon excitations exceeding 1 /zm, which shows 
that ballistic transport of magnetic excitations over rel- 



atively long distances in such systems is indeed possible. 
Various mechanisms to create pure (no charge currents) 
spin currents have been proposed, including spin pump- 
ingi^ and the creation of spin voltage by means of the spin 
Seebeck effect.*^ It has also been shown that, using the 
spin Hall effect, it is possible to convert an electric signal 
in a metal into a spin wave, which can then be trans- 
mitted into an insulator We mention here also the fact 
that this type of ballistic transport in magnetic systems 
is not the only kind of nearly dissipationless spin trans- 
port. Other systems that have received wide attention 
recently are topological insulators jJiS where in particular 
edge states in a quantum spin Hall insulator have been 
predicted to be dissipationlesS;'^^ For an overview of other 
possibilities, see a recent review by Sonin.— 

Specifically, in this work we investigate rectification ef- 
fects in nonitinerant one-dimensional quantum spin sys- 
tems. Our motivation to do so lies therein that such 
effects form the basis for a device crucial to spin- and 
electronics, the transistor. In charge transport, rectifica- 
tion effects in mesoscopic systems have received consid- 
erable attention in recent years, both in one-dimensional 
(lD)2i"— and higher dimensional systems.—"— However, 
these studies have so far been limited to the itinerant 
electronic case. Here, we consider insulating spin chains 
(without itinerant charge carriers) and study the analog 
of such rectification effects in the transport of magneti- 
zation; see Fig. [TJ 

The magnetic systems we consider are assumed to have 
an anisotropy A in the z direction in the exchange inter- 
action between neighboring spins in the spin chain. To- 
gether with an offset magnetic field applied to the entire 
system, also in the z direction, this anisotropy will allow 
us to design systems displaying a nonzero rectification ef- 
fect. We consider both ferro- and antiferromagnetic sys- 
tems, using respectively a spin-wave and a Luttinger liq- 
uid description. For both types of magnetic order we find 
that an enhanced exchange interaction in the z direction 
gives rise to the opening of a gap in the excitation spec- 
trum of the spin chain. We can use the offset magnetic 
field to tune the system to just below the lower edge of the 
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gap, so that we get the asymmetry in the spin current be- 
tween positive and negative magnetic-field gradients Ai? 
required to achieve a nonzero rectification current. We 
find here that larger values of the exchange interaction J 
or the exchange anisotropy A lead to a larger gap, and 
hence require stronger magnetic fields. 

For the case of antiferromagnetic order we also dis- 
cuss the case of suppressed exchange interaction in the 
z direction. As it turns out, in this case we need an 
extra ingredient to achieve a sizable rectification effect. 
Here the extra is a "site impurity," a local change in the 
exchange anisotropy, which models the substitution of 
a single atom in the spin chain. We find that, in the 
presence of the offset magnetic field, the leading term in 
the perturbation caused by such an impurity fiows - in a 
renormalization-group sense- to strong coupling for low 
energies. We can therefore, again in combination with 
the offset magnetic field, use it to achieve the rectifica- 
tion effect. In the regime where our calculations are valid, 
we find that the rectification current is quadratic in the 
strength of the impurity, and behaves as a negative power 
of ( A_B / J) . The dependence of the rectification current 
on the anisotropy A is more complicated, and described 
in detail in Sec. IVlIIl 

This work is organized as follows: In Sec. |TT]we intro- 
duce our model, in Sec. IIIII we discuss the rectification 
effect for ferromagnetic systems using the spin- wave for- 
malism. In Sec. IIVI we map the antiferromagnetic (AF) 
Heisenberg Hamiltonian on the Luttinger liquid model 
and describe the different perturbations resulting from 
the anisotropic exchange anisotropy in the presence of 
the offset magnetic field. We then continue in Sec. |V] 
with a renormalization-group analysis to find the most 
relevant perturbations. In Sec. IVII we focus on the case 
of enhanced exchange anisotropy and find the resulting 
rectification effect. Section fVIII contains the analysis for 
a suppressed exchange anisotropy in combination with 
the site impurity, and in Sec. IVIIII we give numerical 
estimates of the rectification effect for the latter system. 



II. SYSTEM AND MODEL 
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FIG. 1. (Color online) (a) Schematic view of the nonitinerant 
spin system. The ID spin chain is adiabatically connected 
(see text) to two spin reservoirs, (b) A magnetic-field gradi- 
ent A_B gives rise to transport of magnetization through the 
spin chain, e.g. via spin waves. Depending on the exchange 
interaction Ai in the spin chain [see Eq. ([l])], part of the spin 
current can be reflected. In a rectifying system, the magni- 
tude of the spin current through the system depends on the 
sign of AB. 



to be constant, we can distinguish two different classes 
of ground states in the absence of an external magnetic 
field, depending on the sign of the exchange coupling 
Jij . For constant < the spin chain has a ferromag- 
netic ground state and we can describe its low-energy 
properties within the spin-wave formalism. For constant 
Jij > the ground state is antiferromagnetic, and, in 
principle, we can determine the full excitation spectrum 
using Bethe-ansatz methodsi^ However, given the com- 
plexity of the resulting solution it is convenient to use a 
different approach and describe the system and its low- 
lying excitations using inhomogeneous Luttinger liquid 
theory, which is what we will do in this work. 

We will study two different scenarios for the exchange 
interaction in the spin chain. In the first scenario 



The system we consider consists of a one-dimensional 
spin chain of length L, adiabatically connected (see be- 
low) to two three-dimensional spin reservoirs. This sys- 
tem is described by the Heisenberg Hamiltonian 

(ij) i 

(1) 

Here (ij) denotes summation over nearest neighbors, and 
we count each bond between nearest neighbors only once. 
Sf is the a component of the spin operator at position 
fi. Jij denotes the exchange coupling between the two 
nearest neighbors at fi and fj. The non- negative vari- 
able Aij is the anisotropy in the exchange coupling. Re- 
garding the ground state of the spin chain, assuming Aij 



A > 1. 



(2) 



In general, such a constant anisotropy in the spin chain 
will open a gap in the excitation spectrum, which we can 
use to design a system that displays a nonzero rectifica- 
tion effect. Furthermore, it will turn out that when the 
constant A satisfies < A < 1 we need a spatially vary- 
ing Aij to achieve a sizable rectification effect. Therefore 
in the second scenario we will consider a site impurity at 
position jQ. This perturbation, in which two adjacent 
bonds are altered, describes the replacement of a single 
atom in the spin chain. It can be modeled as 

A„- = A -f A'S,,j-i [S^^,, + (5,,,o+i) and A < 1. (3) 

Here Si is the Kronecker delta function. 
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A the excitation spectrum of the magnons, the bosonic 
excitations of the system, has the fohowing form 



FIG. 2. (Color online) Schematic view of the system under 
consideration (bottom) and the band structure for the ferro- 
magnetic exchange interaction (top). The colored areas on 
the left and right in the top picture depict the filled states in 
the reservoirs. The position of the bottom of the band in the 
left reservoir depends on Ai?. The striped area in the middle 
denotes the band of allowed states in the spin chain. It is seen 
that spin transport through the spin chain can be large for 
AS < 0, but only small for AB > 0. 



The last term in the Hamiltonian ([T]) is the Zeeman 
coupling caused by an external, possibly spatially vary- 
ing, magnetic field Bi. This term defines the z axis as the 
quantization axis of the spin. In all cases that we will con- 
sider the total magnetic field Bi can be split into a con- 
stant part and a spatially varying part, Bi — Bq + ABi. 
Here the constant offset magnetic field i?o > is applied 
to the spin chain and the two reservoirs. ABi is constant 
and equal to —AB in the left reservoir and goes to zero in 
the contact region between the left reservoir and the spin 
chain. We now define a rectifying system as a system in 
which the spin current Ig satisfies Is{AB) ^ ~Is{~AB); 
see also Fig. [T] To quantify the amount of rectification, 
we use the rectification current Ir{AB), defined as 



Ir{AB) 



h{AB) + h{-AB) 



(4) 



III. RECTIFICATION IN FERROMAGNETIC 
SPIN CHAINS: SPIN WAVE FORMALISM 

To illustrate the mechanism behind the rectification 
effect we first consider a ferromagnetic spin chain with 
constant anisotropic exchange interaction in the z di- 
rection [Jij = J < and Ay = A > 1 in Eq. ([T])]. 
The spin chain is adiabatically connected (see below) to 
two three-dimensional ferromagnetic reservoirs Rl and 
Rr (see Fig. [1]) and we assume initially that only the 
the constant magnetic field Bq is present. For temper- 
atures T <^ gfiBo/kB the density of spin excitations is 
low, which allows us to use the Holstein-Primakoff trans- 
formation to map the Heisenberg Hamiltonian on a set 
of independent harmonic oscillators In the presence of 
a constant magnetic field Bq and a nonzero anisotropy 



fiLOk = g^iBBo + 2 (A - 1) \J\s + \J\s{ka)^ 



(5) 



with a the distance between nearest-neighbor spins and 
s the magnitude of a single spin. The magnons have 
magnetic moment —g^sz. The magnon distribution is 
given by the Bose- Einstein distribution UBi^^k)- In order 
for the spin chain and reservoirs to be adiabatically con- 
nected, the length of the transition region Lt between 
spin chain and reservoir must be much larger than the 
typical wavelength of the excitations. From Eq. ([5]) we 
can now see that for ferromagnets this requirement be- 
comes Lt 3> 27r-\/ J s / g^BB^a (we ignore the anisotropy 
in this estimate). For Js = IO/cb J, Bo = 0.5 T the 
requirement amounts to Lt 3> 25a. 

From Eq. ([5]) it follows that magnons in the spin 
chain require a minimum energy of gfiBBo + 2 (A — 1) Js, 
whereas the magnons in the reservoirs (assuming that 
A = 1 in the reservoirs) have a minimum energy of 
gfiBBo- The effect of applying a magnetic field —AB to 
the left reservoir is to create the nonequilibrium situation 
in which the distribution of magnons in the left reservoir 
is shifted by an amount — y/isAB. From Fig. [2]it is seen 
that shifting the magnon spectrum up {AB < 0) allows 
the magnons from R^ to be transported through the spin 
chain. Because of the asymmetry of the distributions in 
Rl and Rr a large (negative) spin current will flow in 
this situation. When we shift the spectrum in R^ down 
{AB > 0), the spin current is blocked by the gap in the 
excitation spectrum of the spin chain, hence only a small 
(positive) spin current will flow. To determine the rectifi- 
cation current we use the Landauer-Biittiker approach^ 
to calculate the spin current through the chain. The total 
spin current is given by Ig =^ Il^r. — Ir^l, where 

Il^r{AB) = -gfiB / —nB{LOk)v{k)T{k), (6) 

and Ir^l is defined analogously. The AB-dependence 
will be shown to be in the limits of integration. Here v{k) 
is the group velocity of the magnons, v{k) = dujk/dk, 
and T{k) is the transmission coefficient through the 
spin chain. For this system the transmission coefficient 
T{k) = 1 as long as the magnons are not blocked by the 
gap in the excitation spectrum of the spin chain. In the 
absence of AB the magnon spectrum in the reservoirs 
reaches the upper edge of the gap in the spectrum at 
wave vector fcp = •\/2(A — l)/a'^. We can incorporate 
the shift in the magnon spectrum in the left reservoir 
by changing the lower boundary in the integral for the 
spin current to fcmin,L — max{0,fcL}, where we defined 

kL — Re \/k^ + ^AB and ^ = g^B/{Jsa?). For tem- 
peratures T such that T ^ sJ/kB, so that we can set 
the upper boundary to infinity, we then have the limits 
of integration (max {0, fc^} , c») for the current from Rj^ 
to Rr. For Ir^l we have the limits (maxj/co, ^i?} , oo). 
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FIG. 3. (Color online) Spin current as function of applied 
magnetic field for the parameters J = lOks 3, s — 1, g = 2, 
Bo = O.f T, T = 100 mK, and a = 10 nm. The spin cur- 
rents for different anisotropics all saturate at the same posi- 
tive value; this is the point where all the magnons incoming 
from the left reservoir are transmitted. This maximum cur- 
rent is on the order of 10^ magnons per second. Since the 
typical magnon velocity v = Ja^ko/h is on the order of 10^ m 
s~^, and assuming length L « 1 fim for the spin chain, this 
is indeed within the single magnon regime. 



scribes the spin chain as the fermionic Hamiltonian: 

i i 

+ (4ic,+i - 1/2) (clc - 1/2 



Ho + Hb + Hz 



(8) 



We initially assume that the magnetic field satisfies 
g^isB <^ J and that A <C 1, so that we can use pertur- 
bation theory to describe Hb + Hz ■ We will ultimately 
want a description of the system by its bosonic action. 

Considering first Hq, and restricting the Hamiltonian 
to low-energy excitations, we can take the continuum 
limit and linearize the excitation spectrum around the 
Fermi wave vector kp = 7r/(2a), where a is again the lat- 
tice spacing, to arrive at an effective (1 -I- l)-dimensional 
field theory involving left- and right-moving fermionic ex- 
citations. To this purpose we replace c\ — a}^'^^\x), 
X^i a — / dx and Aj — >■ A(x). Here x — ia. After intro- 
ducing left- and right-moving fermions, -0^(2;) — ip\^{x) + 
we carry out a bosonization procedure'^^''^^ using 
the following operators: 



where kn = Re ^AB] . The resulting spin current is 
then 



ii{ko,kR} 



dk 



2ak 



(7) 

where a = sa^J. The spin current has been plotted in 
Fig. [3] for realistic material parameters. 



IV. RECTIFICATION IN 
ANTIFERROMAGNETIC SPIN CHAINS: 
LUTTINGER LIQUID FORMALISM 



We now consider the case of an antiferromagetic spin 
chain consisting of spins- 1/2 and adiabatically connected 
(for the antiferromagnetic system this means Lt ^ 
2iTkp^ = 4a, see below) to two three-dimensional anti- 
ferromagnetic reservoirs. We will show that it is possible 
to map both the spin chain and the two reservoirs on the 
Luttinger liquid model. 

We start with the description of the spin chain and 
come back to the description of the reservoirs at the end 
of this section. To model the spin chain'^- we use that in 
one dimension we can apply the Jordan- Wigner transfor- 
mation to map the spin operators onto fermionic opera- 
tors: 5*+ ^ c|e'''^5;-~'=]^^ and 5f ^ c\c., - 1/2. This 
allows us to rewrite the part of Hamiltonian ([TJ that de- 



(9) 



Here r = L,R, = Tlj a-nd fc^ is the Fermi wave vector 
for r- moving particles (see below). (I){x) is the bosonic 
field related to the density fluctuations in the system as 
dx4>{x) — —■K[pii{x)+pL(x)\. Here, 0{x) is its conju- 
gate field, satisfying [(f){x),dx'0{x')\ = i-K5{x — x'). We 
have left out the Klein operators in the fermionic creation 
operators, since they cancel in the subsequent perturba- 
tion theory. Using the aforementioned operations we can 
transform the Hamiltonian into the quadratic action 



27ri^ 



dV 



u 



.(10) 



Here and in the following we use the notation r — (x, t) . 
For the free action u = Ja/h and K — \, these param- 
eters will change due to the presence of the Hb and Hz 
terms. 

The effect of the inclusion of the magnetic-field term 
Hb is twofold: In the absence of any backscattering and 
umklapp scattering the field will introduce different den- 
sities of left- and right-moving excitations in the spin 
chain. The effect of this is to change the Fermi wave 
vectors in Eq. ^ for the respective particles to fc/? — 
TT/i2a) + J^gpB {Bo - AB) and fc^ = 7r/(2a) + ^gMsSo- 
This does not affect the bosonized form of Hq, but will 
have an effect on the bosonized form of Hz- Furthermore, 
the spatial dependence of the magnetic field makes that 
the spin chain is now described by S[(l>] = 5*0 [(^] + Sb[4>], 
where: 



Sb^ = 



5MB 



dV</)(f)a^B(r). 



(11) 
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We use the specific form of the magnetic field 

d^B{f) AB6{x - L/2), (12) 
which corresponds to the magnetic field described in Sec. 

m 

In the next step we derive the bosonic representation 
of the Hz termi^i^ In the continuum limit the z compo- 
nent of the spin operator is given by the normal ordered 
expression 

+2PI{x)M^) + ^l{x)ijR{x) : . (13) 
The interaction term contains terms that transfer ap- 
proximately 0,2kp or Akp momentum. Around half fill- 
ing and for a constant exchange anisotropy, e.g., the case 
in Eq. ([5]), conservation of momentum requires that the 
only terms that can survive are the ones that transfer ap- 
proximately or 4fci;^ momentum, the 2kF terms are sup- 
pressed by rapidly oscillating ex:p{±2ikFx) terms. The 
terms that transfer small momentum give rise to terms 
proportional to {dx(j>{r))'^ , and hence change the param- 
eters u, K in Eq. ([TU)) . It is not possible to determine the 
exact values from the linearized theory, they can, how- 
ever, be determined from the Bethe-ansatz solutioui^ 
Important here is that for A ^ 1 we have K ^ 1/2. The 
Akp term (umklapp scattering term) becomes 
aJ 



Svs [4>] 



(2™) 



d'^rA(a;) cos [A(j){f) - Apx] , (14) 



where p = ^gp^B [Bq ~ AB/2). We have neglected a 
constant term cx pa inside the cosine. 

If a spatially varying anisotropy is present, as in Eq. 
dS]), the 2/cj?-backscattering terms do not necessarily van- 
ish in regions where A (a;) varies on a scale of order 
a. The bosonization of the 2kF terms then requires 
some care. If we naively use the continuum form of 
Eq. ([5]), this term contains infinities after the bosoniza- 
tion. Therefore it has to be normal ordered, "^^ which 
we can do using Wick's theorem and the contraction 

Tprix)ipl,{x + a) = —r ^2Trai ^r,r' ■ Sincc the SfSfj^i term 
contains four-fermion operators the normal ordering will 
yield not only additional constants, but also two fermion 
operators, as can, for instance, be seen from the typical 
term ^^^^^{x)^^ r{x)iI)^j^{x + a)ipL{x + a), which becomes 

: '(/;l?,(^)V'-R(a;)V'l?.(2; + «)V''l(x + a) : 

+ : ip''nix)ipL{x + a) : ^^(a;)^'"^^ -|- a). (15) 

The latter term is equal to —^jif^ ■ ^Ri^)'^Lix + 
a) :. There are four such terms; together they give the 
backscattering term 



aim 



(27ra)2 

X {sin [2(l){r} - 2px] + sin [2(l){r) - 2p{x + c^\}) 



(27ra)2 ^ 

X {sin [20(r) — 2px] + sin [20(r) — 2p{x 
The completely normal-ordered term becomes 

SbsA4>] = tP^ I dV(-l)-/'^A(x) [ad^iT^f 



The latter can be seen by expanding the normal ordered 
term around x and invoking the Pauli exclusion princi- 
ple. From simple dimension counting it is clear that the 
dimension of this operator is 2 -I- i^. Hence it is always 
irrelevant, and we will neglect it in the following. The 
total action describing the spin chain is then given by 

= S-o [4>] + Sb + 5us [0] + ^Bs [0] . (18) 

Now we can distinguish between the two scenarios in 
Eqs. ^ and For the constant anisotropy in Eq. ^ 
the ^Bsi^] and S'bsn[0] terms vanish because they are 
proportional to the rapidly oscillating (— 1)^/"^. Hence 
only the bulk umklapp scattering term, given by 



Ai 



(27ra)2 



dxdt cos [40(f) - -ipx] , (19) 



remains. Here Ai = aJA. The spin chain in Eq. ([2]) 
is thus described by the action S'o[(^] -I- S'busM + '5's[(?!>]. 
As we will discuss in the next section, the bulk umklapp 
scattering term renormalizes as 2 — AK, hence it is irrel- 
evant for if > 1/2 (A < 1) and relevant for if < 1/2 
(A > 1). If this term is relevant, as it is for the parame- 
ters in Eq. ^ , it opens a gap in the excitation spectrum 
of the spin chain, which, as we will show in Sec. I VII can 
be used to achieve the rectification effect in a similar way 
as for the ferromagnetic system in Sec. IIIII To wit, if we 
tune Bq such that it lies just below the upper edge of the 
gap, for /S.B > there can be no spin transport, since 
there are no states available for transport in the chain, 
whereas for Ai? < the states above the edge of the gap 
are accessible, and transport is made possible. 

Next we discuss the case with spatially varying ex- 
change anisotropy, Eq. We start out by noting that, 
since the bulk umklapp scattering term is irrelevant for 
this scenario, the spin chain is not in a gapped state in 
equilibrium. The effect of the applied magnetic field is 
then to move the spin chain away from half filling. As 
we will show later, in the current scenario this doping is 
required in order to achieve a nonzero rectification cur- 
rent. In the case of Eq. ([3]) the backscattering term 
vanishes everywhere, except in the region where A(x) it- 
self varies on a short lengthscale. For the specific form of 
the anisotropy of Eq. ([3]) , the action resulting from the 
site impurity is 



5- 



a2JA' 



{2nY 
sin 



di9j;{sin 24i{xQ,t) - 2p{xo - -) 



2^{xo,t)-2pixo + -) }, 



(20) 



where, here and elsewhere, dxf{xo,t) should be read as 
dxfix,t)\x=xo- Adding to this term the local umklapp 
scattering term caused by the site impurity, the total 
action coming from this impurity becomes 



6 



S"! = y I - [A^ COS 4(/)(xo ,0 + ^2 sin 4(/«(a;o , t)] + cr [Ag cos 2(j){xo , + A^ sin 20(a;o , t)] 
+a [A45:r(/)(xo,t) cos2(?!)(a;o,t) + \\dx(j){xQ,t) sin24>{xQ,t)\ 



(21) 



which contains terms caused by umklapp scattering (pro- 
portional to A2'^), terms that may be called offset- 
induced backscattering terms (proportional to Ag'*"), and 
terms that describe combined forward- and backscatter- 
ing (proportional to A4'^). In the expression for the action 
S'i[0] we have defined a = pa = ^gfis (Bq - AB/2) = 

' gtiB(Bo-AB/2) 



K 



Here we have identified the UV cut- 



ofi' of the theory ujc with u/a. Furthermore, the prefactors 
are given by 



A! 
A| 



2A 
2A 



= A(-l)^o/'^ 
= A(-l)^"/" 



cos4p (xq - f ) + cos4p (xo -I- f )] , 
sin4p (xo - f ) + sin4p (xq -h f )] , 
cos 2p [xq — f ) + cos 2p (xq 
sin2p (xo — f ) + sin 2(0 [xo - 



-!)]- 
§)]• 

a, ^^^^ 

Where A = aJA'. Furthermore, A4' — —A3' . In 
this scenario the spin chain is thus described by So[(j>] + 

Finally we return to the description of the spin reser- 
voirs. As we show in detail in Appendix \^ we can 
describe the low-energy excitations of the reservoirs by 
the quadratic Luttinger liquid action, Eq. (ITUl) . The 



effective Luttinger liquid parameters Ur, Kr of the three- 
dimensional reservoirs can be determined by mapping its 
dynamic susceptibility onto that of a Luttinger liquid, 
using the nonlinear sigma model, resulting in 



ViJa/h if^ = 7r/(4V3). 



(23) 



This means that we can describe the reservoirs by letting 
M — > u(a;) and K — > K{x) in Eq. ([T0|) . where 



u{x), K[x) 



u,K for a: e (-L/2,L/2), 
Ur.Kr for X i (-L/2,L/2). 



(24) 



V. RENORMALIZATION GROUP 
TREATMENT 

We start the analysis of the antiferrromagnetic spin 
chain by studying the scaling behavior of the spin chain, 
allowing us to determine which perturbations will be 
most relevant in the low-energy sector. We perform 
the renormalization-group (RG) calculation in momen- 
tum space,-^ assuming there is a hard natural momen- 
tum cutoff Aq in the system. In the RG procedure the 
cutoff A(Z) = Aqc"' is decreased from A(Z) to A(/ -I- dl). 
For the RG procedure we consider the partition function 
written in terms of the action in imaginary time. As is 
customary we split the field ^(r) contained in this ac- 
tion up in a fast and a slow part, </'(r) = (jy'{r) + 4)'^{r), 



where the fast part contains Fourier components with 
momentum between A(^ + d^) and A(Z), and the slow 
part contains the components with momentum less than 
A(Z -t- dl). The RG procedure then consists of integrat- 
ing out the fast modes, and subsequently restoring the 
original cutoff in the action, in order to find an effective 
low-energy action with the same couplings, but different 
coupling constants. This allows us to find the relevant 
(increasing in magnitude under a decrease of the cutoff) 
and irrelevant (decreasing in magnitude under a decrease 
of the cutoff) couplings. For completeness, we mention 
the renormalization equations for the constant anisotropy 
A (see Ref. [Ml) 



dK _ 8\l 1 
"dT " ~ (27rft)2 ^ 

du 8XlK^ 1 



KC 



K, 



(25) 



dl {2TTfh? u*^"' 

dp 

d/ 

f .(2-4^)A.. 

We have omitted several terms Aa, which are a number 
of order 1. The different constants are given by [here 
f = A(a;,MT)"^ is dimensionless] 

Ck^\J d'rcos [Apx] fVo(f)e-«^^l(""^ 

Cu^-\j d^f{x^ - u^f^) cos [Apx] Jo(r)e-8^^i(^"). 

(26) 

Here Jo(f) is the zeroth order Bessel function and 
Fi{f) = /g dg-!— Both integrals converge and 
are of order 1. From the last RG equation it follows 
that for K < 1/2, which as we have seen before corre- 
sponds to A > 1, the perturbation caused by the bulk 
umklapp scattering grows under renormalization. This 
corresponds to the opening of a gap in the excitation 
spectrum of the spin chain. The magnitude M of this 
gap has been calculated analytically using Bethe-ansatz 
methodSj-"^- and is given by 

M ^ TT sinh ^ 1 

J V ^ cosh[(2n+ I)7r2/2H' ^ ' 

n— — oo '-^ / ! 1 

where v — acoshA. For A > 1 this gap is exponentially 



small, M w 47rJexp 
A — cx) it becomes M 



-7rV((2[2(A-l)]i/^^ 
=s J[A-2]. 



and for 
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If we now assume to be in the regime where K > 1/2 
and add the site impurity described by 5'i[0], we get the 
following additional equations: 



dA^/" 
dl 


-(1- 


AK) A^/' 


~ -1- 44 '^l 


~ ^33 f-^/S: 


dxl"' 

dl 




K)\l"' - 


i 23 




dxl'" 

dl 


= (1- 


K)\l''' - 


i 24 





(28) 

where we have defined the second-order 



terms Tnm 



AJ^Am + A^Am and A„{; 



^r/Aj'i^'' sin4pxo + A^^"^ cos4/9a;o^ Ai. Here 77 = ±1. 
The different constants used here are given by 



ZTT'^hu / r 



K 

■n'^hu 
K 



dfJo(f)e-3^^^(^), 

d2fJo(f)e-2^^i(^), (29) 
dV cos [4(p/A)x] Jo(f)e-3^^i(^), 
d^f^sin [4(p/A)x] Jo(r)e-3^^i(^). 



Again, all integrals converge and are of order 1. From 
the form of the equations it follows that we can ignore 
contributions that are of second order in the couplings, 
and we see that the most relevant couplings are the A3''' 
and A4' terms, which to first order scale as 1 — iiT. For 
e (f/2, 1), the regime in which we operate in the case 
of Eq. (jSj, these terms thus grow in magnitude under 
a decrease of the cutoff. Because of the extra dx(t>{xQ^t) 
proportionality in the A4''' terms, one would expect the 
effect of these terms on the spin current to be smaller 
than that of the A3''' terms. However, since the A3''' 
terms have an additional a in front of them, in the regime 
Bq « AB/2 these are suppressed, so that they could be- 
come comparable to the A^'*" terms. Therefore we will 
consider the effect of both types of coupling when calcu- 
lating the spin current for the system with anisotropy in 
the exchange interaction as in Eq. ([3]) in Sec. IVHl 



VI. ENHANCED ANISOTROPY 

As we have shown in Sec. |Vl the enhanced exchange 
anisotropy of Eq. ([2]) gives rise to a bulk umklapp scat- 
tering term in the action of the spin chain, and the chain 
is described by the sine-Gordon modelj^ This model has 
been well studied, and is known to give rise to two pos- 
sible phases depending on the value of the chemical po- 
tential g^sB: For values of the chemical potential lower 



than the gap M the system is a Mott insulator; when the 
chemical potential is increased to values above this gap, 
the system undergoes a commensurate-incommensurate 
(C-IC) transition and becomes conducting. 

To simplify the sine-Gordon model we replace (/)m('^) = 
2(f){r) and, in order to keep the commutation relation be- 
tween the two fields unchanged, 0M{r) = 9{r)/2. The 
umklapp scattering term then reduces to a backscatter- 
ing term (a two-particle operator), and the effective Lut- 
tinger parameter becomes AK. The fermionic form of the 
resulting new action is known as the massive Thirring 
model. At the Luther- Emery point, K ~ 1/4, (Ref. 
[ssh the new fermions whose action is given by the mas- 
sive Thirring model are noninteracting, and we can di- 
agonalize the quadratic part of the action by a Bogoli- 
ubov transformation, which gives rise to two separate 
bands of fermionic excitations with dispersion e±(fc) = 
±yJ{ukY + M'^. UK ^ 1/4 there are residual four- 
fermion interactions present. However, it can be shown^ 
that, independent of the initial interactions, near the C- 
IC transition the strength of these interactions vanishes 
faster than the density of the fermions, so that the inter- 
actions become negligible. The new free fermions are not 
the original fermions; instead they correspond to soli- 
tonic excitations of the original action. These solitons 
have fractional magnetic moment — g/iB/2. We finally 
can relinearize the excitation spectrum around the point 
M to arrive at a new Luttinger liquid in terms of the 
fields ipAii'T), SAiir) with parameter Km = AK. 

To calculate the dc spin current through the system we 
use^ that in linear-response theory the dc spin current is 
given by Is = GAB, where the conductance G is given 
by 



G = 



-I hm 



(30) 



and G^^{x,x' ,iOn) is the time-ordered Green's function 
in imaginary time 



G^4x,x',uj^)^ / dre"^"^ {Tr(j){x,T)cl){x',0)}s„. 
Jo 

(31) 

Here a;„ is the Matsubara frequency. At zero tempera- 
ture, and assuming infinitesimal dissipation in the reser- 
voirs, the w — >■ limit of this Green's function can be 
determined for the inhomogeneous system including the 
two reservoirs?^ and is given by 2^^. The effect of the 
entire mapping of the original sine-Gordon model onto 
the new free Luttinger liquid can be captured here by 
replacing gfiB — >■ 9fJ'B/'2 and Kr AKr, so that the 

conductance in the conducting phase is G = Kr ^^^^ . 
Following Ref. |4ll we conclude that excitations that are 
injected at energies above the Mott gap are transported 
through the chain, giving rise to the aforementioned con- 
ductance, whereas excitations injected at chemical po- 
tential lower than the gap are fully reflected. Assuming 
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that giisBo ~ M, this gives rise to the spin current 

Is{AB) = Kr^^^^ABOi^AB), (32) 

where Q{—AB) is the unit step function. Since spin 
transport is absent for AB > 0, we have Ij.{AB) — 
Is{AB)/2. In these calculations we assumed that the 
length of the spin chain L — >■ oo, so that we can neglect 
tunneling of solitons, and we have assumed zero temper- 
ature. Finite size and temperatures are known to give 
corrections to the conductance.— 



Here, Io{AB) is the spin current through the systems in 
the absence of Si [(j)] , given by the well-known expression 



I,iAB)=K^-^AB, 



(35) 



and I^s{AB) describes (negative) contributions to the 
spin current due to Si[(j)], which we will derive next. 

The contribution to the spin current resulting from 
a (T^p- cos[20(O, t)] term (which describes offset-induced 
backscattering) is given by 

g/iscr^ (A? r)^ 1 
haiAB) = \ ''^^ KrAoiAB), (36) 



VII. SUPPRESSED ANISOTROPY 



where 



To determine the rectification current resulting from 
the exchange anisotropy as given in Eq. ([3]) we need to 
calculate the spin current in the system given the action 
Sy)] = So[(p] + Si[(j)] + 5'b[(/)]. For simplicity we assume 
here that the impurity is located at xq = 0. We ignore 
the bulk umklapp scattering term, since we have shown 
in Sec. |V] that the contribution of this term is irrele- 
vant for the parameters used here. From the RG analy- 
sis it also followed that the most important terms are the 
offset-induced backscattering terms and that in regions 
where Bq « AB/2 the effect of the terms describing com- 
bined forward- and backscattering and the offset-induced 
backscattering terms can become comparable, due to the 
extra a in front of the latter terms. Therefore we need to 
calculate the spin current due to both types of contribu- 
tion. We will show that the rectification effect appears 
in the contributions to the spin current that are second 
order in the coupling constants. 

We calculate the spin current using the Keldysh tech- 
nique,^'* We assume that at t = —oo the system is de- 
scribed by the action S'o[0] + S'b[0], and that the per- 
turbation Si[(j)\ is turned on adiabatically. From con- 
servation of spin it follows that the expression for the 
spin current in the Luttinger liquid is given by /s('') — 
— ^^dt4>{r)- The spin current can then be calculated 



/, = - 



TT \ 6J(f) 



Ao{AB) 



K4: 



r(l + 2K) 



IR \1r\ 



-2+2K 



-2|7hI 



(37) 



Here we introduce the dinicnsionless parameters Ag ^ = 
Xi^/{hhJc) and — K^gnBAB / {hwc). As before, de- 
notes the UV cutoff of our theory, given by Uc = u/a. 
In Sec. |V] it was determined that the backscattering 
term scales as 1 — K under renormalization. To im- 
prove our result we should therefore not use the bare 
coupling Ag^, but instead the renormalized coupling. 
Since we assumed an infinitely long chain, and consider 
zero temperature, the renormalization group procedure 
is stopped on the energy scale determined by the mag- 
netic field, gfisAB. We can account for this by replacing 

AR -1 + ^ 

^'s.R ^ ^^hu ^t,R- P°i"* ^^'^^ deter- 

mined the spin current resulting from the backscattering 
term. By repeating the previous calculation with both 
the Ag- and the Ag-proportional terms included it is eas- 
ily seen that the (Ag) -proportional contribution to the 
spin current is also given by Eq. ([36|) . with (Ag^)^ re- 
placed by (Ag ^)^. The two cross terms proportional to 
AgAg do not contribute to the spin current, since they 
cancel each other. 

The calculation of the contribution to the spin cur- 
rent due to a term ^^9a;</'(0, t) cos 20(0, t) that describes 
combined forward- and backscattering proceeds along the 
same lines; we refer again to Appendix IB] for the details. 
The resulting contribution to the spin current is given by 



Here {(f>^{r))g is the average of the field (f>{r) over the 
Keldysh contour with respect to the action 5" [(/)], where 
the ± denotes that the field is located on the positive, 
respectively negative, branch of the contour. The right- 
hand side contains the functional derivative of the parti- 
tion function of the system with respect to the generating 
functional J(r), which is given in Eq. (jB2l) . The details 
of the calculation are given in Appendix |B] here we sum- 
marize the results. We find that there are two fundamen- 
tally different contributions to the spin current. 



IsiAB) ^ loiAB) + iBsiAB) 



ha{AB) = 



KrAi{AB), (38) 



where 



Ai{AB) 



T{2 + 2K) 



(39) 



Again, according to the RG analysis, we have to replace 
the bare coupling with its renormalized value, in this 



case, A^^^ 



gfj-B^B 



-l+K 



XI Like with the backscat- 



(34) tering terms, we can easily determine the effect of the Af 
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and A| terms combined, 
to (A4) is given by Eq. 



The spin current proportional 
with (A4 ^) replaced by 



(At 



and the two cross terms cancel each other. 



Finally, we need to consider the cross terms between 
the A3 and A4 terms, such as, for instance, A3A4. Us- 
ing the results from Appendix [O it is easily seen that 
these always vanish, since they are all proportional 

±2iU"(0,t)-'^'''(0,t')] 



to < rA0"(O,i)e 
backscattered current is then given by 

5Ms K, 



/bs(AB) 



where A,;,eff = 



1/2 



0. The total 

,ff^i(Ai3)}, 
(40) 



Equation (PH)) is the main result of this section. From 
the explicit form of Ao{AB) and Ai{AB) [see Eqs. (071) 
and (|39|) ] it is clear that the impurity flows to strong 
coupling for low AB, as was expected from the RG anal- 
ysis. We note that, since both Aq{AB) and Ai{AB) are 
odd in AB, one could naively think that the impurity 
does not contribute to the rectification current, which 
requires the spin current to have a component that is 
even in AB. However, from the explicit form of the A^^off 
[see Eq. ([^^ for the bare couplings] it follows that these 
couplings also contain parts that are proportional to AB. 
Furthermore, the part of the spin current that is propor- 
tional to Aq{AB) is proportional to tr^, which also con- 
tains a AB. Physically, these terms are caused by the 
fact that an incoming excitation sees a slightly different 
impurity depending on the energy with which it comes 
in. The spin current Eq. (j40[) therefore has components 
even in AB, which contribute to the rectification current. 
We are now also in a position to explain why it is needed 
to move the system away from half filling in order to ob- 
tain a nonzero rectification current. If we set Bq — in 
Eq. (|22p. it turns out that the A^^eff's and a contain only 
terms that are even in AB, so that the spin current is 
odd in AB. In contrast, when Bq ^ 0, we create terms 
in the Ai^eff's and a that are odd in AB, thereby causing 
a nonzero rectification current. 



VIII. ESTIMATE OF RECTIFICATION 
CURRENTS 

In this section we will give numerical estimates of the 
rectification current for realistic experimental parame- 
ters. Experimental results^^ show that the exchange cou- 
pling in certain spin chains can be on the order of J w 10^ 
K. The effect of a different exchange interaction is illus- 
trated in Fig. m In the analysis of these results one 
must keep in mind that, since our perturbative calcula- 
tion of Ibs{^B) diverges for AB — >• 0, it breaks down for 
field differences |AB| < |AB*|, where AB* is the field 
such that max[/Bs(AB*),/Bs(-AB*)] = /o(AB*). In- 
stead of showing the apparent divergent behavior, the 
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FIG. 4. (Color Online) Normalized rectification current 
Jr/|-^o| (main figure) and backscattering current /bs/I-^^oI (in- 
set) as a function of the applied magnetic-field difference AS 
for different values of the exchange interaction J. Param- 
eters are K = 0.63, Bo = 750 mT, and A' = 0.5. The 
black dots in the main figure denote the values AB* for 
which max[|JBs(AB*)!,|JBs(-AB*)|] = |/o(AB*)|. As ex- 
plained in the text, our perturbative results are not valid for 
\AB\ < |AB*1. 



backscattering current must go to zero for AB < AB*. 
With this in mind. Fig. 2] shows that a smaller exchange 
interaction J gives rise to a larger rectification current 
at equal AB. The maximum value of Ir{AB) will also 
be reached at a higher value of AB. In order to get the 
largest possible rectification current it is thus required 
to use a material with an exchange interaction as small 
as possible, with the constriction that it must be large 
enough to yield its maximum at reasonable values of AB. 

In Fig. [5] we show the dependence of the backscat- 
tering and rectification current on the Luttinger liquid 
parameter K. The behavior for smaller K is similar to 
the behavior shown for smaller J: The maximum rectifi- 
cation current is increased, but occurs at a higher value 
of AB. Indeed, we know that both /^(AS) and Ir{AB) 
obey a power-law dependence of AB with negative expo- 
nent. Since the modulus of this exponent decreases for 
increasing K, the behavior is as expected. 

Figure [5] shows the dependence of the backscatter- 
ing current on the applied magnetic field Bq. A larger 
amount of doping clearly corresponds to larger possible 
rectification current, again at the price of a higher re- 
quired magnetic field AB. This can be understood by 
realizing that Bq determines to a large extent how much 
of the impurity the incoming excitations see at low ener- 
gies AB Bq, as follows from the a dependence in the 
Ag'*" terms in the action Eq. (|21l) . 

In Ref. H it has been shown that N parallel uncoupled 
AF spin chains connected between two three-dimensional 
AF spin reservoirs, each one carrying a spin current 
Is{AB), give rise to an electric field 



7V^M^(O,cos20,-sin20) 



(41) 



Here the spin chains are assumed to extend in the x di- 
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FIG. 5. (Color Online) /bs/|-^o| as a function of the applied 
magnetic field difference AB for different values of K. Pa- 
rameters are Bo = 750 mT, J = 100 K, and A' = 0.5. See 
the caption of Fig. |4]for the explanation of the black dots. 



is by spin accumulation. By applying again an ac driv- 
ing field to the left reservoir it is possible to measure an 
accumulation of spin in the right reservoir, since trans- 
port is asymmetric with respect to the sign of AS. We 
consider again 10** parallel spin chains with K = 0.6, 
consider i?o = 750 mT, and an amplitude AS — 43 mT 
for the driving field. For /S.B{t) < the spin current is 
always zero. For AB{t) > there is a nonzero spin cur- 
rent. Assuming that the spin current is 10% of the un- 
perturbed spin current [the value at AB{t) — AB] over 
the entire range {0,AB), the rate of spin accumulation 
is about 10~^^ JT~^s~^ « 10"'^^ magnons per second. We 
also note that, contrary to the electric case, where charge 
repulsion prevents a large charge accumulation, there is 
no strong mechanism that prevents spin accumulation in 
this scenario. 
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FIG. 6. (Color online) /bs/I-'^oI as a function of the applied 
magnetic- field difference AB for different values of Bo. Pa- 
rameters are K = 0.63, J = 100 K, and A' = 0.5. See the 
caption of Fig. |4]for the explanation of the black dots. 



rection, and z is the quantization axis as before. Also, 
— y/r, and cosip = z/r. We can 



+ z' 



sin( 



use this electric field to measure the rectification current 
Ir{AB). To illustrate the method we assume K = 0.6, 
Bo = 750 mT, and TV = 10"*. We apply the time- 
dependent field AB{t) ~ ABcos{u]t) to the left reser- 
voir. From Eq. (^(11) it follows that, if we trust our per- 
turbative calculation of Ibs up to the value AB* « —43 
mT, where |/o(AB*)| fa |/bs(AB*)|, the difference in 
magnitude between \Is{AB*)\ and \Is{-AB*)\ is on the 
order of 10% of the unperturbed current |/o(Ai3*)|. From 
Eq. (HIT) , and assuming r = 10~^, it then follows that 
the difference in voltage drop between two points (0, r, 0) 
and (0,0, r) is Vs ~ 10~^^ V, which is within experi- 
mental reach. Here we note that as long as the driving 
frequency w <C J/fi, our calculation of the spin current 
in the dc limit remains valid. Another constraint is given 
by the spin- relaxation time r^, which is typically on the 
order of 10~^ s."*^ This allows us to use frequencies in the 
MHz- GHz range. 

Another possibility to observe the rectification effect 



IX. CONCLUSIONS 

In this work we have proposed various realizations of 
the rectification effect in spin chains, consisting of ei- 
ther ferromagnetically or antiferromagnetically coupled 
spins, adiabatically connected to two spin reservoirs. For 
both ferromagnetic and antiferromagnetic spin chains we 
find that the two crucial ingredients to achieve a nonzero 
rectification current are an anisotropy in the exchange 
interaction in combination with an offset magnetic field, 
both in the z direction. An enhanced anisotropy creates 
a gap in the excitation spectrum of the spin chain. Us- 
ing the magnetic field to tune the chemical potential of 
the excitations to just below the lower edge of the gap 
then allows us to achieve a large rectification current for 
realistic values of the applied magnetic field. For antifer- 
romagnetic coupling and suppressed exchange interaction 
in the z direction we find that a uniform anisotropy is not 
sufficient to achieve a sizable rectification effect; instead 
we use a spatially varying anisotropy, which we attain in 
the form of a site impurity. Away from half filling this 
impurity fiows to strong coupling under renormalization, 
which allows us to achieve a sizable rectification effect for 
realistic values of the applied magnetic field. We have 
also proposed several ways to observe this rectification 
effect. 
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Appendix A: Determining Luttinger Liquid 
parameters from the nonlinear sigma model. 

In this section we calculate the magnetic susceptibil- 
ity of a three-dimensional antiferromagnet, denoted by 
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Xq/sC?)- In the continuum limit, the magnetic suscepti- 
bility is defined as 

X^?(9) = (^)' ;^ / d'^e^'-^ {TrS^{s)Spm , (Al) 

where s = (t, r)-^, q = (a;„, fc)"^, and q ■ s = uJuT — k ■ r. 

We start from the isotropic Heisenberg model in three 
dimensions, described by Eq. ^ with A^- = 1 and Bi — 
0. We introduce the unit vector field nf = Si/S. This 
vector can be written as a staggered Neel vector part 
Hi plus a fluctuating part li, both of which are slowly 
varying on the scale a. 



nl = (-lyn^^^ I - aH^ + ah « {-lyn. + ak. (A2) 

In the continuum limit it can be shown^* that the imag- 
inary time action describing this system can be approxi- 
mated as S[n, I] ~ S{)\n, l\ + S'/[n, Z], where 

3 



S'o[n, /] 



h 

2^ 



dS 



c ^ 



Si[n,l] = I2ha / dS 



/ + i 



1 



{n X drfi) 



,(A3) 



.(A4) 



Here, a = J/{8ah), c — \/3Ja/h, and 7 = 4\/3a'^- From 
a substitution of Eq. (IA2|) in Eq. (|A1[) it is clear that we 
can write 



3D/ N /.9AiB\^ 1 /■ 1. 



(r,Z,(s)/^(0))5, (A5) 



where the Green's function is calculated with respect 
to the action S[n,l]. The Green's function is deter- 
mined from the partition function Z ~ J DnDlS{n^ — 



1) exp y^fi ^S[n, l]j . We can perform the functional in- 
tegration over I in Eq. (jASP to arrive at 



1 

he 



V 7 

X /(n X a^n)^ (s) (n x (0) 



So 



where we replaced r — > cr. Here, the average is with re- 
spect to Sq [ft, I] . Assuming the magnetization is ordered 
along the z axis, we write ft = {ai,a2,TV) . The a^'s 
describe thermal fluctuations perpendicular to the direc- 
tion of order, and H = Vl — ct^. Also, a = {ai,a2)^. 
For low temperatures, we can then rewrite the action of 
Eq. dMI as 



27 



d^^ 



EC 

1=1 



(A7) 



From this and Eq. (|A6p it then follows that the magnetic 
susceptibility in the gapless direction is given by 



icq)' 



(A8) 



To compare this relation to the magnetic susceptibility 
in one dimension, we use X^^{q) = o,^X^^{q)- Since we 
know^ that in the Luttinger liquid model 



(A9) 



we have 



VSJa/h Kr = 7r/(4V3) (AlO) 



as effective Luttinger liquid parameters of the three- 
dimensional reservoirs. 



Appendix B: Calculation of the spin current in the 
Keldysh formalism 

In this appendix we will show how to calculate the 
spin current /^(Ai?) in a system described by the action 
S[(j)] = So[(j)] + Si[(t)] + SB[(t)]- According to Eq. the 
spin current is given by 



»,=± 



^1 



We start by deriving the contribution to the current due 
to a cos 2(/){0, t) term. In this case the partition func- 
tion can be written in vectorized form as 

Z[J{r)] = I Dcfcxp - i y" dVdV 

f) ■ Go(f, r')${f') - 2iS{r^ r')${r) ■ QJ{r') 

(B2) 



X exp 



where we have defined 
J{r) 

Q 

Go(r',r') 



dt ^ 77C0S [20''(O,i)] 
-00 ^=± 



J(f)/2 
G++(r,r') G+-{r,n 



(B3) 



Here Gq'^ is the contour-ordered Green's function 



So 



(B4) 



We will also need the retarded, advanced, and Keldysh 
Green's functions G^'^^^^ {f,r'), which are related to 
the contour-ordered Green's functions in the usual way. 
Performing the linear transformation (j)' {r) — <j){r) — 
i / dr^Go(r, r')(5"^ J(r') shows that the partition function 
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can be factorized into two parts, the first one containing 
the quadratic part of the action, the second containing 
the backscattering term. The spin current is therefore 
made up of two contributions. The first, Iq, comes from 
the quadratic action, and describes the spin current in 
the absence of any backscattering. The second, /3a, orig- 
inates from the backscattering term. Using the exphcit 
form for the magnetic field Eq. ((T^ it turns out that the 
expression for the unperturbed spin current is given by 



The contour-ordered correlation function can be deter- 
mined from the inhomogeneous Luttinger liquid model^ 
and the fluctuation-dissipation theorem, and can be writ- 
ten in terms of f as 



T,[^+{0,T)-q}-{0,0)Y) =i^lnHf + l]. (BIO) 



dt'G^{f;~L/2,t'). (B5) 



Since this is exactly the expression for the spin cur- 
rent in the linear-response regime,— which yields Iq — 



r h 



AB, we can infer that 



With this expression we can rewrite Eq. (jB9[) as a func- 
tion of 7fl' and K only 

MAB) = -J.(^^7«l7«r^^'"e-I-I. (Bll) 

Next we discuss the contribution of the 

^d^(j}{0,t)cos2<j){0,t) term from Eq. dH]). We 

again perform the linear transformation (f>{r) — )■ 4''{r), 
which yields 



dt'iG^[f--Ll2,t') = 



irKrt 



(B6) 



Idn. — 



^^^a, / di'<(f;0,O 



To calculate the contribution from the backscattering 
term we perform the functional derivative of the partition 
function with respect to J{r). The resulting expression 
is 



'3a 



x(sin[2(/.'+(0,t')-27<'])5, (B7) 
where 7 = Kr ^^^^'^^ . Expanding this term in Af gives 



X {d^4'+{0,t')sm [20'+(O,t') - 27i'])g(B12) 

There is also an additional term in the expression for 
the spin current introduced by the linear transformation, 
which is proportional to dx J dtG§{0,t;r"). Equation 
(jB6p . however, shows that this term is zero. We can 
expand Eq. (|BT2)) in A^^, 

2gfiB{Xtfa^ 



-dt / dt'iG^{f;0,t')dt" 

<J —OQ 



ha = i- 



-dt 



dt'iG{l{f;0,t')dt" X 



7] (sin [2(j)'+{0,t') - 2jt'] cos [20"'(O,i") - 27t"])^^ correlation function. 



X J2 V{d.'cl>'+{0,t')sm [2cf,'+{0,t') - 2^t'] 

xdx"(f>"'{0, t") COS [2(f>"^{0, t") - 27i"] .(B13) 
We now use the results from Appendix [Cl to evaluate the 



Since contour ordering reduces to time ordering for two 
fields on the same branch of the Keldysh contour, and 
because the system in equilibrium is time-reversal invari- 
ant, the term with the fields at times t' and t" on the 
same branch averages to zero. Switching the variable 
of integration to t = t' — t" and using Eq. jBH} it is 
seen that the contribution from the backscattering term 
is given by 

guBcr'^ (AS r)^ 1 
/3a = — ^-T^ KrAoiAB). (B8) 

TT* as 

Here — a/u can be interpreted as the cutoff in the time 
domain. As before, we identify the cutoff in time with 
the cutoff in frequency: as = i^jT^, so that we can write 
the expression for the backscattering current in terms 
of the dimensionless parameters A3 ^ = X^/{fiuJc), Jr — 
KrgfiB^B /(hujc), and f — r/as- Here 

Ao{AB) = - / dfsin[27flf]e ^ ^'^ ' ■ ' " ' • / so 
^ J -00 

(B9) 



OA (T,0+(O, (0,i')) 



1 + if 



^0 (uas)2 \1 -I- r 



(B14) 

We also use that {Tc(j)+{0,t)(j)-{0,t'))g^ = 0, so that 
the spin current is given by 

TT* as 

where the dependence on the parameters of the spin chain 
and the applied magnetic field difference is contained in 
the function Ai{AB) 

X {dxC^'{Q, r)5,(/.'(0, 0))s^ sin [27;^f] , (B16) 

where we again used the fact that terms containing cor- 
relators with both fields on the same branch of the con- 
tour vanish. Straightforward manipulations show that 
Ai{AB) is given by 

A, [AB) = IB. iT^r e-^l-^^l ■ (B17) 
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Appendix C: Correlators 

In this appendix we show how to calculate a correlation 
function of the form 



T,nAM^,,U)e'^^^^^^^'-''A . (CI) 

/ So 



First of all, we can pull the derivatives out- 
side of the expectation value, so we can calculate 

Tcn,(/.(xj,tOe'^^^'^(="^'*^^) instead and take the 

derivatives afterward. For a lighter notation we prove 
the required results in discrete space, making the con- 
nection = (jiiri). We use the main result of Gaussian 
integration^ 



(C2) 



Given the partition function 

Z -- 



we can then calculate, for instance, {<t>i4>j) So letting 
9? w U work on both sides of Eq. (IC2I) . To calculate 

Jnjm Ij— ^ ^ 



the required correlator, Eq. (|C1|) . we replace j ^ j + b 
in Eq. ()C2p . where b consists of the components bi = 
J2j ^j^ij- By letting 9j„|j^q work on both sides of this 
equation we arrive at 



(C4) 



Here we sum over double indices. We also used the fact 
that is a symmetric matrix, e.g., = G~^, which 
is true since it contains contour ordered correlation func- 
tions. By letting d'f U work on Eq. (|C2[) we get 



[<PnC^^e' = [G-i + G-,]kG;,lk] o-^", 

(C5) 

Both equations above are valid when J^j = 0; other- 
(C3) wise the correlators vanish. 
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